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Abstract 



We have searched for the pah production of first generation scalar leptoquarks 
in the eejj channel using the full data set (123 pb^^) collected with the D0 
detector at the Fermilab Tevatron during 1992-1996. We observe no candi- 
dates with an expected background of approximately 0.4 events. Comparing 
the experimental 95% confidence level upper limit to theoretical calculations 



' of the cross section with the assumption of a 100% branching fraction to eg, 

■ we set a lower limit on the mass of a first generation scalar leptoquark of 225 

O ■ GeV/c^. The results of this analysis rule out the interpretation of the excess 



of high events at HERA as leptoquarks which decay exclusively to eq. 
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Leptoquarks (LQ) are hypothesized exotic color-triplet bosons which couple to both quarks and 
leptons. They appear in extended gauge theories and composite models ^ and have attributes of 
both quarks and leptons such as color, fractional electric charge, and lepton and baryon quantum 
numbers. Leptoquarks with universal couplings to all flavors would give rise to flavor-changing 
neutral currents and are severely constrained by studies of low energy phenomena Q. Therefore, 
only leptoquarks which couple within a single generation are considered here. The HI and ZEUS 
experiments at HERA have reported an excess of events at high in e~^p collisions One 
possible interpretation of these events is resonant production of first generation leptoquarks |Q] at a 
mass near 200 GeV/c^. To date, no excess has been observed in e^p collisions Q. A straightforward 
leptoquark explanation then requires the leptoquarks to decay to eq with a branching fraction of 
100% 1]. 

The CDF and D0 collaborations have published the results of searches for first generation 
leptoquarks in pp collisions using data collected prior to 1994 [^j. The HI and ZEUS experiments 
at HERA have published lower limits on the mass of a first generation leptoquark which depend on 
the unknown LQ-e-q coupling ||5[0|. Experiments at LEP have searched for leptoquarks in decays 
of the Z boson @. 

This Letter describes a search for the pair production of first generation scalar leptoquarks in 
the eejj + X final state using 123ib7 pb~^ of data collected at the Fermilab Tevatron with ^/s = 1.8 
TeV during 1992-1996. The D0 detector and data acquisition system are described in detail in 
Ref. The detector consisted of three major subsystems: a uranium/liquid-argon calorimeter; 
central tracking detectors, including a transition radiation detector; and a muon spectrometer. 

A base data sample of 101 events with two electrons and two or more jets was selected. Electrons 
were identified by their longitudinal and transverse shower profiles in the calorimeter and by the 
fraction of their energy deposited in its electromagnetic section. The electrons were required to 
be isolated from other energy depositions and to have pseudorapidity \r]\ < 1.1 or 1.5 < \r]\ < 2.5. 
In addition, at least one electron was required to have a matching track in the central tracking 
detectors and to satisfy ionization requirements in the tracking chambers and transition radiation 
detector. Jets were reconstructed using a cone algorithm of radius TZ = \J (A0)^ + (A?])^ = 0.7, 
where is the azimuthal angle, and required to have |?7| < 2.5. The electrons were required to 
be separated from jets by 7^ > 0.7. Electrons were required to have transverse energy E?p > 
20 GeV and jets to have Ej. > 15 GeV. The kinematic quantities were calculated using the vertex 
determined by the electrons. Events whose ee invariant mass lies between 82 < Mee < 100 GeV/c^ 
{Z boson region) were rejected. The efficiency of the electromagnetic trigger used to collect the 
base data sample exceeded 99% for the leptoquark mass range addressed by this analysis. 

Monte Carlo (MC) signal samples were generated for leptoquark masses between 120 and 260 
GeV/c^ using the ISAJET |1C| event generator and a detector simulation based on the GEANT |11] 
program. Leptoquark production cross sections were taken from the recently available next-to- 
leading order (NLO) calculations of Ref. [^]. The primary backgrounds to the eejj decay mode 
are Drell-Yan + 2 jets production (DY), ti production, and misidentified multijet events. Monte 
Carlo samples for the DY events were generated using ISAJET. The DY cross section normalization 
was fixed by comparing the MC events with Z + 2 jets data in the Z boson region. Top quark 
events were generated using the HERWIG |l^ program at a top quark mass of 170 GeV/c^ and 
all dilepton final states were included. The D0 measured it production cross section of 5.5 it 1.8 
pb at a top quark mass of 173.3 GeV/c^ was used |14|. The multijet background was estimated 
from a sample of events with four or more jets in which the probability for two jets or photons 
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FIG. 1. St distributions for background (solid line histogram), data (solid circles), and 200 
GeV/c^ leptoquark MC events (open triangles). 



to be misidentified as electrons was weighted by the number of jets in the event which passed the 
electron Et and r] requirements. This misidentification probability was calculated from a sample 
of events with three jets to be (3.50 it 0.45) x lO"'^ for an electron with a reconstructed track and 
(1.25ib0.13) X 10~^ for an electron without a reconstructed track. The errors on these probabilities 
reflect a slight dependence on the jet Et and rj. The signal and background samples were required 
to pass the kinematic selection criteria that defined the base sample. 

To search for leptoquarks, a random grid search method ||T^ was used to optimize cuts on the 
data and MC samples. Consistent results were obtained using a neural network [16|. Two criteria 
were used to optimize event selection, one designed for discovery and the other for limit setting. 
Once it became clear that the data did not contain evidence for leptoquark production, the limit 
setting criterion of a maximum number of signal events for a fixed number of background events 
was adopted. The background level chosen was 0.4 events, corresponding to a 67% probability that 
no such events would be observed. 

The set of cuts which optimally separates signal from background was determined using a 
systematic search over a grid of possible cuts with the choice of grid points determined by the 
distributions of signal MC events. Many sets of selection criteria were explored including com- 
binations of kinematic quantities and mass-related variables, such as transverse energy and the 
two ej invariant masses. A cut on a single, relatively simple variable, St = + Hj., where 
Hj^ = Eq} + Ej? and Hij, = X^jets -^T' satisified the limit setting criterion. Approximately 0.4 
background events are expected for St > 350 GeV. No events remain in the base data sample after 
this St cut is applied. Figure |l| shows the St distribution for the base data sample, the predicted 
background, and a MC sample of 200 GeV/c^ LQ events. The highest value of St seen in the data 
is 312 GeV. 

For the neural network (NN) analysis, a three layer feed- forward network was constructed 
with two inputs, Hj, and Hj,, and one output, the NN discriminant -Dmv- Figure |2|(a) shows the 
expected distribution in Hj^ vs. Hj, for the 200 GeV/c^ MC signal sample; Figs. |2|(b) and (c) show 
the same distributions for the predicted background and the base data sample. The network was 
trained using the 200 GeV/c^ LQ sample and the background samples described above. -Dmv has 
a range between (background) and 1 (signal). Figures §(a-c) show contours corresponding to 
three values of -Dmv- A background of 0.4 events is obtained by requiring Datat > 0.95. After 
application of this cut, no events remain in the data. As the NN and 5^ analyses provide nearly 
identical sensitivity, only the simpler St analysis was used for the cross section limit described in 
this Letter. 

The background was estimated for St > 350 GeV and is given in Table | for the three sources. 
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FIG. 2. Hj^ vs. Hj, for (a) 200 GeV/c^ LQ events, (b) predicted background, and (c) base 
data sample. The curved lines correspond to Dmn = 0.5, 0.8, and 0.95 (from left to right). The 
area of a box is proportional to the number of events in the bin, with the total number of events 
normalized to 123 pb^^. 



TABLE I. Background contributions from individual sources. 



Background Source 


Number of Events 


DY 


0.18 ±0.04 


tt 


0.11 ±0.04 


Multijet Misidentification 


0.16 ±0.02 


Total 


0.44 ± 0.06 




The total estimated background is 0.44 ± 0.06 events where the error includes both statistical and 
systematic uncertainties. Included in the systematic error are the uncertainties in the jet energy 
scale, particle identification efficiency, ti production cross section, and luminosity, and the effects 
of the choice of parton distribution function and renormalization and factorization scale ;U, gluon 
radiation, and MC statistics. 

Modeling of the St distribution for high mass DY events was checked by studying and 
Hj, separately using data and MC events in the Z boson mass region. The average value of Hj~ 
for high mass DY events (which provide most of the DY background) is approximately 250 GeV, 
corresponding to an Hj, of approximately 100 GeV for St = 350 GeV. The distribution of Hj, for 
high mass DY events is expected to be similar to that of Z ± 2 jets events. Figure ^ shows the H?p 
distributions for Z + 2 jets MC and data. In the region corresponding to the St cut for high mass 
DY events {H^ ~ 100 GeV), the agreement is good. 

To investigate the background further, constrained mass fits were performed on the events in 
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FIG. 4. 5't' vs. 3C fit mass distributions for (a) background, (b) 200 GeV/c^ leptoquarks, and 

(c) the base data sample. Tlie area of a box is proportional to the number of events in the bin. 

(d) Mass distribution of the events in the base data sample (solid circles), expected background 
(solid line histogram), and 200 GeV/c^ leptoquarks (hatched histogram). The inset plot shows 
these distributions for events with St > 250 GeV. 



the base data sample, on background samples, and on the 200 GeV/c^ leptoquark signal sample. 
The 3C mass fit was based on the SQUAW [jl^ kinematic mass fitting program and required the two 
ej masses to be identical. Use of the fitting program improves the mass resolution by approximately 
10% over a simple calculation of the ej invariant masses. Figures Q(a-c) show St as a function 
of the fit mass for the estimated background, 200 GeV/c^ leptoquark events, and the base data 
sample. The distribution from the data agrees with that of the background. The two highest mass 
events have low values of St and so are unlikely to be leptoquark events. Figure ^(d) shows the one 
dimensional mass distributions for the same samples. Inset in Fig. ^(d) are the distributions after 
a cut on St > 250 GeV. As can be seen, the data are consistent with the background prediction. 

The dielectron identification efficiency was determined to be (73 it 4)% using a sample of 
Z ee+2 jets events. The overall signal detection efficiency is 9-37% for leptoquark masses of 120— 
250 GeV/c^ (Table II). We set a 95% confidence level (CL) upper limit on the cross section a using 
a Bayesian approach with a flat prior distribution for the signal cross section. The statistical and 
systematic uncertainties on the efficiency, the integrated luminosity, and the background estimation 
were included in the limit calculation with Gaussian prior distributions. The resulting upper limit 
on the cross section is shown in Fig. |^ together with the NLO calculation of Ref. |12|, and the 
results are listed in Table ^ The intersection of our limit curve with the lower edge of the theory 
band (/i = 2Mlq) is at o" = 0.068 pb, leading to a lower limit on the mass of a first generation 
scalar leptoquark of 225 GeV/c^. For a branching fraction of 100% to eq, this is our lower limit; 
inclusion of additional channels will provide increased sensitivity in the case where leptoquarks 
also decay to uq. 
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TABLE II. The signal detection efficiency, the 95% CL upper Umit on the production cross 
section, and the fi = 2Mlq NLO cross sections from Ref. |12|. 



Leptoquark Signal 95% CL NLO Theory 

Mass Efficiency Upper Limit Cross Section 

(GeV/c^) {%) (pb^ (pb) 

120 8.7 ± 1.4 0.170 3.8 

160 20.8 ±3.0 0.113 0.68 

200 31.2 ±3.8 0.078 0.16 

225 35.7 ±4.3 0.068 0.068 

250 37.2 ±4.5 0.066 0.030 




%0 190 200 210 220 230 240 250 
M^Q (GeV/c^) 

FIG. 5. Upper limit on the leptoquark pair production cross section for 100% decay to eq. 
Also shown is the NLO calculation of Ref. [jl2| where the central solid line corresponds to /i = Mlq, 
and the lower and upper dashed lines to /i = 2Mlq and /x = Mlq/2, respectively. 



In conclusion, we have excluded the interpretation of the HERA high events as first gener- 
ation scalar leptoquarks which decay exclusively to eq, as expected in chiral models with no extra 
fermions or intergenerational mixing. Using the NLO cross section calculation of Ref. [|l^ with 
/i = 2Mlq, and assuming a 100% branching fraction to eq, the 95% CL lower limit on the mass of 
a first generation scalar leptoquark is 225 GeV/c^. 
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and T.G. Rizzo for many useful discussions. We thank the staffs at Fermilab and collaborat- 
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(Prance), State Committee for Science and Technology and Ministry for Atomic Energy (Russia), 
CNPq (Brazil), Departments of Atomic Energy and Science and Education (India), Colciencias 
(Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and CONICET 
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